Testing for double inflation with WMAP 
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With the WMAP data we can now begin to test realistic models of inflation involving multiple 
scalar fields. These naturally lead to correlated adiabatic and isocurvature (entropy) perturbations 
with a running spectral index. We present the first full (9 parameter) likelihood analysis of double 
inflation with WMAP data and find that despite the extra freedom, supersymmetric hybrid poten- 
tials are strongly constrained with less than 7% correlated isocurvature component allowed when 
standard priors are imposed on the cosomological parameters. As a result we also find that Akaike 
& Bayesian model selection criteria rather strongly prefer single-field inflation, just as equivalent 
analysis prefers a cosmological constant over dynamical dark energy in the late universe. It appears 
that simplicity is the best guide to our universe. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Our universe shows evidence of complexity and, at 
the same time, great simplicity. Our universe appears 
entirely consistent with being a "double-de Sitter sand- 
wich" - radiation and matter dominated phases caught 
between two de Sitter phases at low and high energies 
respectively. 

Recent work [J, has shown that a cosmological con- 
stant provides a better fit than dynamical dark energy 
to current CMB and SNIa data if one computes the 
Bayesian evidence or uses information criteria for model 
selection. In this paper we will show that, at least within 
a class of double, hybrid inflation models, the same is 
true for the early universe. One might envisage various 
infrared-ultraviolet dualities to explain such behaviour. 

Despite this apparent "asymptotic blandness " there 
is interesting tentative evidence to the contrary. The 
WMAP data show unusual characteristics such as "oscil- 
lations" which may disappear with more data or may 
be the first signs of new physics. Similarly there is evi- 
dence for a feature in the power spectrum £| which can 
easily be produced by the subtle dynamics of multiple 
light scalar fields during inflation. 

Multiple light fields during inflation automatically 
widens the narrow predictions of single-field inflation 
for now there are multiple entropy perturbations 
which are, in general, correlated to some degree with the 
standard adiabatic mode |l q -[27 | . Correlations are pro- 
duced when the valley of the effective potential is curved 
[l8| and this also leads to non-gaussianity [2|| . Since the 
effective masses of the various fields typically depend on 
the vacuum expectation values of the other (dynamical) 
fields these are time-dependent and can cause violations 
of standard slow-roll conditions and spectral indices for 
the perturbations which run with scale plj . 

This is a crucial aspect of this present work because 



previous analyses of correlated adiabatic and entropy 
(isocurvature) perturbations have always assumed power- 
law spectra for all the perturbations I2 I. IMlM Wj \ . 
When applied to the WMAP data they found that with 
standard priors on cosmological parameters the degree of 
correlation allowed is small (although see H^j). Allow- 
ing running of the spectral index, at least in the super- 
symmetric hybrid models we study, does not change this 
conclusion. 

The code we have developed allows us to numerically 
study any inflationary model without approximation (ex- 
cept in the treatment of spinodal/tachyonic instabilities) 
and builds on that used in [2lJ • Future work will consider 
more fields where sharp features can occur, something 
which does not occur in the two-field double-inflation 
models we study here. 



II. GENERAL FORMALISM 

We consider two minimally coupled scalar fields, 4> and 
X, with an effective potential V(4>i X)- Our main inter- 
est is the case of double inflation in which two stages 
of inflation are realized. General scalar metric per- 
turbations about the flat Friemann-Lemaitre- Robertson- 
Walker background can be written as (see e.g. [lql 



ds 2 = -(1 + 2A)dt 2 + 2a(t)B^dx l dt 
+a 2 (t)[(l-2%h)S, 



(1) 



2E Aj \dx l dx J 



where a(t) is the scale factor. The comoving curvature 
perturbation in our two-field system is then given by 



H(tft5<f> + x Sx) 



x 



(2) 



where H = a/ a is a Hubble rate, and 8<j) and S\ are the 
perturbations of the fields <fi and x> respectively. 
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The perturbation equations are given in Refs. MM 
and one can numerically evaluate the power spectrum, 
Vn = (k 3 /2n 2 )\K\ 2 , at the end of inflation j^H (here k is 
comoving momentum). In the multi-field system we also 
need to account for the spectra of isocurvature pertur- 
bations, Vs, and correlated adiabatic and isocurvature 
perturbations, Vc (see Ref. for their definitions). 

The quantity re defined by rc = Vc / VV-rVs is the mea- 
sure of the strength between adiabatic and isocurvature 
perturbations. 

The system possesses several model parameters associ- 
ated with the potential. We assume slow-roll conditions 
apply \4>\ <C |3ff</>| and \x\ <C |3i?xL for the initial con- 
ditions of background fields, so that the 4> and x terms 
are neglected. Then the initial conditions of <j) and \ are 
determined by 0i n and Xxa (the subscript "in" denotes 
the initial values). We perform the likelihood analysis 
over the initial conditions Xin- 

Note that the number of inflationary model parameters 
depends explicitly on the inflaton potential and typically 
requires at least three parameters in the context of double 
inflation. 

We impose the condition that the total number of e- 
folds during inflation must exceed Nt = 50 to solve flat- 
ness and horizon problems. We find the cosmologically 
relevant perturbation modes with comoving wavenum- 
bers k and numerically evolve the background and all 
perturbation equations through inflation, giving us the 
three power spectra (V-r.,Vs,Vc), as described in |2l|. 

It is important to solve the perturbation equations 
without approximation right up to the end of inflation, 
since the curvature perturbation is not necessarily con- 
served after Hubble radius crossing Q , unlike the case of 
single-field inflation. 

The resulting data: Pi = logV (ki) given at a wave 
number of Xi = logfc;, i = 1, ..to, are optimally fitted 
with a polynomial function Pn t (x) — a + bx + ex 2 by 
minimizing x 2 = J2i( a + + cx 1 — Pi) 2 ■ 

For each set of parameters we derive the best-fit coef- 
ficients a,b,c for each of the three power spectra. It is 
worth mentioning that the coefficients b and c are inti- 
mately linked to the spectral index n s and its running 
of scalar perturbations by the relations n s = 6+1 and 
a s = 2c. We check that our fitting method agrees very 
well with numerically obtained power spectra and is suf- 
ficient to accurately capture any running of the spectral 
index over cosmologically relevant scales. 

We assume, as is standard, that the field <fi decays 
to ordinary matter like photons, neutrinos and baryons, 
whereas the field x decays into cold dark matter (CDM) 
[l6l Us! . In this case the mixing between two scalar 
fields is negligible and the CDM isocurvature perturba- 
tions and correlations remain after reheating. Relaxing 
this assumption will introduce extra parameters into the 
analysis. 

The CMB temperature anisotropics are given in gen- 



eral by 

C e = {An) 2 J k 2 dkA 2 (k,T ), (3) 

where Ag(k, tq) is the £-multipole of the A:-th wavenumbcr 
temperature anisotropy at the present time To. For a 
general set of correlated initial conditions, one has 

A(fe,r ) =P^ 2 A n (k,T ) + P 1 s /2 A s (k,T ), (4) 

where (A^(/c,r in )) = (A|(fe,%)) = 1 and 
(A K (k,r m )As(k,T in )) = r C - Then we get 

C e = (4tt) 2 J k 2 dk[P n Al K + P s Al s + 2r c (A e , n A e , s }] 
= Ce^n + Ce^s + 2Ce,c ■ (5) 

It is possible to obtain the three multipole spectra re- 
quired for any general set of initial perturbations using 
the following simple scheme. Let us denote C(Ai 1 A2) 
as the Ce spectrum obtained with completely correlated 
initial conditions with a given adiabatic spectrum A\ 
and given isocurvature spectrum Ai. A typical Boltz- 
mann code can produce only C(Ai,0) (pure adiabatic), 
C(0,A 2 ) (pure isocurvature) or C(Ai,Ai) (completely 
correlated mixture of adiabatic and isocurvature with the 
same initial spectrum). It is not difficult to see that the 
general spectrum is given by: 

C(A 1 ,A 2 ) = C(A 1 ,0) + C(0,A 2 ) + 2[C(A 12 ,A 12 ) 

-C(0,A 12 )-C(A 12 ,0)], (6) 

where in our case Ai 2 = \[P~c , A\ = \fpji and A 2 — 
\/P~s. One needs therefore five evaluations for each com- 
bination of Pji,S,C- We make use of a modified version of 
the CAMB Boltzmann solver [3(j to evaluate the CMB 
power spectrum by this scheme. 

In addition to 4>\ ni Xin and the inflationary potential 
parameters discussed in the next section, we varied 4 
cosmological parameters: flbh 2 , fl c h 2 , r, Ho; namely the 
baryon and cold dark matter density, the reionisation 
optial depth and the Hubble constant today. We assume 
spatial flatness, so fl\ = 1 — f2f, — f2 c . 

It is well-known that the allowed ranges for these pa- 
rameters has a large impact on the acceptable amount 
of correlated isocurvature perturbations [31|. We choose 
fairly standard priors, allowing the above variables to 
vary in the ranges: r S [0,0.3], H e [50,90] with fi;, 
and £l c both varying over the full unit interval, [0, 1]. We 
choose very wide domains for <p- m and Xia and found that 
the results depended very weakly on these boundaries. 

We then use the first year WMAP TT and TE data 
|32| in our analysis to constrain the various parameters. 
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III. A REALISTIC DOUBLE INFLATION 
MODEL AND LIKELIHOOD RESULTS 

Let us consider a fairly realistic multi-field inflation 
model with potential 

F=j(x 2 -^) 2 + ^W + ^ 2 , (7) 

corresponding to the original version of the hybrid infla- 
tion [23j. This is closely linked with those obtained in 
supersymmetric theories 0, HE 0- which generically 
leads to a very strong correlation between the adiabatic 
and isocurvature perturbations due to the presence of a 
tachyonic instability between the two phases of inflation 
|2l| . In this work we concentrate on the supersymmet- 
ric case with g 2 /X = 2. Then we have three potential 
parameters: A, M and to, which are constrained by our 
likelihood analysis. 

We can have two stages of inflation for the potential 
(|7|l depending upon the model parameters. One corre- 
sponds to the stage with <j> > <p c = M/g driven by the 
slow-roll evolution of 4> during which the potential is ap- 
proximately described by V ~ M 4 /4A + m 2 (j> 2 /2. An- 
other inflationary stage is the one with cj) < 4> c driven by 
the field \ with a tachyonic instability. 

When the condition M > to is satisfied, then M 4 /A > 
to 2 2 , and so the Hubble rate is roughly constant with 
a value H = \/2tt / '3AM 2 / m p , around <p — <p c (here m p 
is the Planck mass). We can estimate the condition for 
double inflation by estimating the effective masses of the 
two fields, i.e., ~ to 2 and to 2 ~ g 2 (f> 2 — M 2 . Double 
inflation occurs when both of the masses of the two fields 
are smaller than H, which gives the condition 

M 2 > mm pX / 3A/27T , (8) 
(M/m p ) 2 > 3A/2tt. (9) 

We are mainly interested in the double inflation sce- 
nario in which the second stage of inflation occurs after 
the symmetry breaking. Since m x is smaller than H 
around = C , the field x is hardly suppressed during 
the first stage of inflation, unless tfi is not too much larger 
than (j) c . 

On the other hand, when (M/m p ) 2 <C A, the field x 
is exponentially suppressed for <fi > (j) c and rapidly wa- 
terfalls toward the global minimum of the potential af- 
ter symmetry breaking. This corresponds to the original 
version of the hybrid inflation without a second stage of 
inflation |33j. 

In this case the homogeneous mode of x can be van- 
ishingly small relative to its fluctuations, so the analysis 
using linear perturbations is not fully trustworthy. In 
our work the linear perturbation equations in Ref. |2lj 
are used to evaluate the three power spectra at the end 
of inflation. While the system is stable for the parameter 
range in which double inflation occurs, we found a strong 
numerical instability for perturbations in the tachyonic 



instability region when the field \ is strongly suppressed 
before symmetry breaking. Thus the latter case is effec- 
tively excluded from our analysis. In this case we need 
to account for the effect of diffusion using e.g., a Fokker 
Planck equation |35ll37| . but we do not consider this here. 

In order to constrain the double inflation model given 
by Q), we perform a likelihood analysis over 9 parame- 
ters: 5 inflationary (M, A, to, 0; n , Xin) an d 4 cosmological 
(H , n b h 2 , n c h 2 7 r). 

A grid-based analysis over all 9 parameters would re- 
quire a great deal of time and computing resources, and 
would still lead to very coarse sampling of the param- 
eter space. Instead we conducted the analysis using a 
Markov Chain Monte Carlo (MCMC) approach. We ran 
independent chains on different HPC facilities and used 
the Gelman and Rubin statistic to test for con verg ence 
and mixing of our MCMC chains, as discussed in |32jl38l |. 

Our 2d likelihood plots show two different results for 
the I and 2-cr contours • The filled contours are com- 
puted by binning the MCMC chains, and drawing con- 
tours around points where the likelihood x 2 has dropped 
to 0.32 (l-cr) and 0.05 (2-cr) respectively 1 . 

On the other hand, the unfilled line contours show the 
regions which contain 68% (1-cr) and 95% (2-cr) of all the 
points in our chains (after burn- in phases are removed). 
We define the burn-in point for a chain to be the place 
where the x 2 of the chain drops below the global median 
X 2 for the first time, as in [40| . 



A. Inflationary parameters 

In Fig.^ we show the 2-dimensional likelihood 
plots for various combinations of dimensionless in- 
flationary parameters: (M/m p , (M/m p ) 2 /A) and 
(m/M, (M/m p ) 2 /X). From the left panel it is clear that 
the 2-cr likelihood area is clustered in a small region 
around A <~ (M/m p ) 2 . The square of the effective mass 
of x relative to H 2 is given as |m 2 /i? 2 | ~ 0.5|(c6/</> c ) 2 — 1| 
for A - (M/m p ) 2 . Therefore | m x \ is smaller than H 
for 4> < (j) c , which means that the second stage of 
inflation occurs after the symmetry breaking. When 
A ~ (M/m p ) 2 , the condition JHJ translates into 
M > \/3/2ir m. From the right panel of Fig. ^ one finds 
that m/M varies in the range 0.2 < m/M < 0.7, in 
which case the condition for the first stage of inflation 
is satisfied. Therefore double inflation actually occurs 
within the 2-cr likelihood region of Fig.^ 

In Fig.[21we plot the likelihood constraints for the ini- 
tial values of the scalar fields. These are also constrained 
to lie in a narrow region in the range I.I < 4>i n /4> c < 1.8 
and 2 x 10~ 3 < Xin/xo < 8 x 10" 3 (here xo = M/y/X). 



1 We use x 2 to denote the standard statistical estimate of likeli- 
hood (— ln£) in order to distinguish it from the square of the x 
field. 
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FIG. 1: 2-dimensional likelihood constraints for the double 
inflationary parameters for the potential (|7J. We show la 
and 2a contour bounds from \ 2 data (shaded contours) and 
the contours which contain 68% and 95% of all the points in 
our MC chains (solid lines, see discussion in the text). 



This means that initial values of <f> close to <p = <ft c are 
favoured. Since m x is much larger than H for (f> 3> (j) c , 
the field x is strongly suppressed for the initial condi- 
tions 4>- m /4>c ^S> 1. This corresponds to the case in which 
the perturbations exhibit violent growth in the tachyonic 
region, thus effectively ruled out in our linear analysis. 
The initial value of x affects the number of e-folds dur- 
ing the second stage of inflation (= A^nd)- We obtain 
smaller A^nd for larger Xin/xo- As we find in Fig-El the 
likely values for the number of e-folds is 50 < A^nd ^ 65 
which corresponds to initial conditions Xin/xo of order 

io- 3 -icr 2 . 

It is rather surprising that the likelihood contours of 
A^nd are clustered in the region with cosmologically rel- 
evant scales. In order to obtain this result we did not 
put any prior for the maximum values of the total num- 
ber of e-folds. We found that it is difficult to satisfy 
the conditions of COBE normalization and suppressed 
isocurvature perturbations unless A^nd ranges in the re- 
gion 50 < A^nd ^ 65. This implies that double inflation 
has a rich and complex structure relative to single-field 
inflation. 



B. Power Spectrum 

In this subsection we consider the contribution of 
isocurvature perturbations to the CMB anisotropics. In 
Fig.^jwe plot observational contour bounds for the ampli- 
tude V-r, and the two ratios Ps/Pn, Pc/Pn- The most 
likely value of V-r is around P-r — 2.5 x 10~ 9 , which is 
similar to the case of single- field inflation |24l l41j . 

The contribution of isocurvature perturbations is re- 
quired to be small relative to adiabatic ones to be com- 
patible with CMB anisotropics. As shown in Fig.[S] the 
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FIG. 2: 2-dimensional likelihood constraints for the initial 
conditions of the field, {(f>i n /4>c, Xin/xo}- 




FIG. 3: Marginalised 1-d likelihood of the number of e-folds 
occurring during the 2nd phase of inflation. The solid line is 
based on the number of MCMC points, while the dotted line 
weighs each point based on its \ 2 value, as described in the 
caption to Fig. 1. 



TT spectrum in the isocurvature dominated case does 
not fit with the WMAP data at all. When isocurvature 
perturbations are comparable in magnitude to the adia- 
batic spectrum (labelled "mixed"), the spectrum shows 
significant deviations from the WMAP data on larger 
scales. We found the 2a bounds: Ps/Pn < 0.004 and 
Pc/Pn < 0.07 in order to be consistent with WMAP. 

In Fig.HJwe plot the observational contour bounds on 
A2nd and the spectral index n s . There are some regions 
in which the spectrum of scalar perturbations is blue- 
tilted (n s > 1) with A^nd ^ 55. Since the power spectra 
generated in the first and second stages of inflation are 
blue- and red-tilted respectively [2l| , it is possible to have 
some suppression of power at low multipoles provided 
that the number of e-folds during the second stage of 
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FIG. 4: 2-dimensional likelihood contours for the amplitudes 
of the power spectra, with adiabatic against isocurvature 
(left) and adiabatic against correlated (right). Despite the 
freedom in allowing running of the spectra index the isocur- 
vature component is severely constrained. 
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inflation satisfies 51 < A^nd ^ 55. We show one example 
of the power spectrum in such a case in Fig.0 Although 
strong suppression around £ = 2, 3 is not easily achieved 
unless the spectrum is highly blue-tilted in this region 
(see e.g. |42j), it is intriguing that this double inflation 
scenario provides a possibility to get a better fit on large 
scales. 



IV. DOUBLE INFLATION VERSUS 
SINGLE-FIELD INFLATION 

A natural question is whether the extra complexity 
and fine-tuning involved in double inflation is actually 
preferred by the data over standard single-field inflation. 
This can be addressed by using the Akaike information 
criterion (AIC) and Bayesian Information criterion (BIC) 
|43j. These two criteria are defined as: 

AIC = -2\nC + 2K, (10) 
BIC = -21nC + K1n N p . (11) 

Here C is the maximum value of the likelihood, K is the 
number of parameters and N p = 1348 is the number of 
WMAP data points. The optimal model minimises the 
AIC or BIC. In the limit of large N p , AIC tends to favour 
models with more parameters while BIC more strongly 
penalises them (since the second term diverges in this 
limit) . BIC provides an estimate of the posterior evidence 
of a model assuming no prior information. Hence BIC is a 
useful approximation to a full evidence calculation when 
we have no prior on the set of models (in this case single 
versus double inflation). In this case, we have no strong 
reason a priori to favour double inflation over single field 
inflation so BIC provides sensible approximation to a full 
evidence calculation. 



FIG. 5: The CMB angular power spectra for five different 
cases: 

(i) our best-fit double inflation model, 

(ii) isocurvature dominating over the adiabatic, 

(iii) the isocurvature is comparable to the adiabatic (mixed), 

(iv) the best fit single-field model with potential 111 21 and 

(v) a model with blue-tilted spectrum (n s > 1) on large scales. 
The spectra are significantly different from the standard one 
when the isocurvature is dominant. 
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FIG. 6: 2-dimensional likelihood contours for the spectral in- 
dex n 3 and the number of e-folds during the second stage of 
inflation A^nd- 
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Model -2m £ AIC BIC 

Double inflation 1428.85 1446.85 1493.70 

Single-field 1430.99 1444.99 1480.43 

TABLE I: The best-fit % 2 (= -2 In £) and best Akaike and 
Bayesian Information criteria (AIC and BIC) for single and 
double inflation. Both criteria favour single field inflation. 

Our double inflation model has 5 inflationary parame- 
ters (M, m, A, <pi n , Xin)- We compare this with a single- 
field scenario with potential 

This has 3 inflationary parameters (A, M, Xin)- There are 
also 4 cosmological parameters, common to both models. 
In Table |U we show the best-fit x 2 and the values taken 
by the criteria for the models we have considered. 

We find that the best-fit value of —2 In C in double 
inflation is smaller than in the case of single-field infla- 
tion. However both the AIC and BIC values for double 
inflation are significantly larger than those in the latter 
case, which suggests that single-field inflation is favoured 
relative to double inflation. In addition one could ar- 
gue that single light-field inflation should theoretically 
be preferred a priori since it does not require fine-tuning 
to achieve more than one field to be light relative to the 
Hubble constant. Adding this prior will further favour 
single-field inflation. 

We have only included WMAP data. Evidence for run- 
ning of the spectral index from WMAP and lyman-a data 
|44| would favour double inflation models in which tilt is 
generic |2ll. However evidence for running is currently 
weak (43 and hence should not affect our conclusions 
significantly. Strong evidence for running in future data 
might change the situation however. 

V. CONCLUSIONS 

In this paper we have studied observational constraints 
on double inflation using the WMAP first year data. The 
model we adopted is the supersymmetric hybrid potential 
given in Eq. J7J). The presence of a tachyonic instabil- 
ity region after symmetry breaking leads to the correla- 
tion between adiabatic and isocurvature perturbations, 
which can significantly alter the CMB power spectrum 
compared to the case of adiabatic perturbations alone. 

Comparing with first year WMAP CMB data we found 
that the correlated isocurvature component can be at 
most 7% of the total contribution which is dominated by 
the adiabatic spectrum. 

We carried out likelihood analysis in terms of 5 in- 
flationary parameters and 4 cosmological parameters. 



The likelihood values of inflationary parameters are clus- 
tered in a narrow region around M/m p ~ 5.0 x 10 ~ 8 , 
A - (M/m p ) 2 and m/M ~ 0.5 (see Fig.0. 

In spite of the large number of freedom of model pa- 
rameters relative to single-field inflation, the parameter 
space of double inflation is severely constrained. This 
comes from the fact that it is not so easy to satisfy all con- 
straints including COBE normalization and sufficiently 
suppressed isocurvature perturbations. 

We also found that the number of e-folds in the sec- 
ond stage of inflation are constrained to lie in the range 
50 < A^nd ^ 65. Loss of power on large scales (relevant 
to achieving suppressed CMB low multipoles) is possible 
when the number of e-folds is around 51 < N^nA 55. 

We also compared double inflation with single-field in- 
flation by using the Akaike (AIC) and Bayesian informa- 
tion criteria (BIC). While the minimum value of x 2 in 
double inflation is slightly smaller than in single-field in- 
flation, the information criteria strongly support single- 
field inflation over the supersynmmetric hybrid double 
inflation models we studied. 

Nevertheless we need to caution that the minimum \ 2 
is still larger than the number of data points N p in cur- 
rent observations. We expect that future high-precision 
data such as the Planck satellite will provide more sophis- 
ticated information to distinguish between double infla- 
tion and single-field inflation. 

In this regard it will be interesting to extend our analy- 
sis to include more fields, so that the matter power spec- 
trum can exhibit sharp features, and to allow more real- 
istic treatment of reheating. Both of these will increase 
the number of inflationary parameters (by about 2 or 3 
each) and it is difficult to imagine them producing smaller 
values of the AIC and BIC as a result. 

It seems likely therefore that single field inflation will 
continue to be the scenario to beat. It is intriguing that 
both the early and late universe seem well-described by 
very simple inflationary stages, and perhaps even two 
pure de Sitter phases. Finding a theoretical basis for this 
perplexing high-energy/low-energy duality may become 
a dominant quest in cosmology in the coming years. 



ACKNOWLEDGEMENTS 

We thank Nicola Bartolo, Pedro Ferreira and An- 
drew Liddle for useful discussions. S.T. is grateful to 
Rome observatory, Universities of Sussex, Queen Mary 
and Portsmouth for their warm hospitality during which 
part of this work was done. B.B. is supported by a 
Royal Society-JSPS fellowship. The analysis was carried 
out on the multiprocessor machines Solent and Vela in 
Portsmouth and the UK national cosmology supercom- 
puter (COSMOS) in Cambridge. 



7 



[1] T. D. Saini, J. Weller and S. L. Bridle, Mon. Not. Roy. 

Astron. Soc. 348, 603 (2004). 
[2] B. A. Bassett, P. S. Corasaniti and M. Kunz, arXiv:astro- 

ph/0407364. 

[3] J. Martin and C. Ringeval, Phys. Rev. D 69, 083515 
(2004); P. Hunt and S. Sarkar, arXiv:astro-ph/0408138. 

[4] P. Mukherjee and Y. Wang, Astrophys. J. 599, 1 (2003); 
N. Kogo, M. Matsumiya, M. Sasaki and J. Yokoyama, 
Astrophys. J. 607, 32 (2004); D. Tocchini-Valentini, 
M. Douspis and J. Silk, arXiv:astro-ph/0402583; 
N. Kogo, M. Sasaki and J. Yokoyama, arXiviastro- 
ph/0409052. 

[5] D. Polarski and A. A. Starobinsky, Phys. Rev. D 50, 6123 
(1994). 

[6] A. A. Starobinsky and J. Yokoyama, arXiv:gr- 
qc/9502002. 

[7] J. Garcia-Bellido and D. Wands, Phys. Rev. D 53, 5437 
(1996). 

[8] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 
(1996). 

[9] T. Chiba, N. Sugiyama and J. Yokoyama, Nucl. Phys. B 

530, 304 (1998). 
[10] T. Kanazawa, M. Kawasaki, N. Sugiyama and 

T. Yanagida, Phys. Rev. D 61, 023517 (2000). 
[11] S. Tsujikawa and H. Yajima, Phys. Rev. D 62, 123512 

(2000). 

[12] A. A. Starobinsky, S. Tsujikawa and J. Yokoyama, Nucl. 

Phys. B 610, 383 (2001). 
[13] J. c. Hwang and H. Noh, Class. Quant. Grav. 19, 527 

(2002). 

[14] S. Tsujikawa and B. A. Bassett, Phys. Lett. B 536, 9 
(2002). 

[15] S. Groot Nibbelink and B. J. W. van Tent, Class. Quant. 

Grav. 19, 613 (2002). 
[16] D. Langlois, Phys. Rev. D 59, 123512 (1999); D. Langlois 

and A. Riazuelo, Phys. Rev. D 62, 043504 (2000). 
[17] R. Trotta, A. Riazuelo and R. Durrer, Phys. Rev. Lett. 

87, 231301 (2001). 
[18] C. Gordon, D. Wands, B. A. Bassett and R. Maartens, 

Phys. Rev. D 63, 023506 (2001). 
[19] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 

64, 123504 (2001). 
[20] D. Wands, N. Bartolo, S. Matarrese and A. Riotto, Phys. 

Rev. D 66, 043520 (2002). 
[21] S. Tsujikawa, D. Parkinson and B. A. Bassett, Phys. Rev. 

D 67, 083516 (2003). 
[22] F. Di Marco, F. Finelli and R. Brandenberger, Phys. Rev. 

D 67, 063512 (2003). 
[23] L. Amendola, C. Gordon, D. Wands and M. Sasaki, Phys. 

Rev. Lett. 88, 211302 (2002). 



[24] H. V. Peiris et al, Astrophys. J. Suppl. 148, 213 (2003). 
[25] M. Bucher, K. Moodley, N. Turok, Phys. Rev. D62, 

083508 (2000); M. Bucher, K. Moodley, N. Turok, Phys. 

Rev. D66, 023528 (2002). 
[26] J. Valiviita and V. Muhonen, Phys. Rev. Lett. 91, 131302 

(2003). 

[27] P. Crotty, J. Garcia-Bellido, J. Lesgourgues and A. Ri- 
azuelo, Phys. Rev. Lett. 91, 171301 (2003). 

[28] N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 
65, 103505 (2002); F. Bernardeau and J. P. Uzan, Phys. 
Rev. D 66, 103506 (2002). 

[29] K. Moodley, M. Bucher, J. Dunkley, P. G. Fcrreira and 
C. Skordis, arXiv:astro-ph/0407304. 

[30] A. Lewis and A. Challinor, Phys. Rev. D. 66, 023531 

(2002) . 

[31] M. Bucher, J. Dunkley, P. G. Ferreira, K. Moodley, C. 
Skordis, Phys. Rev. Lett. 93, 081301 (2004) 

[32] L. Verde et al, Astrophys. J. Suppl. 148, 195 (2003). 

[33] A. D. Linde, Phys. Rev. D 49, 748 (1994). 

[34] E. J. Copeland, A. R. Liddle, D. H. Lyth, E. D. Stewart 
and D. Wands, Phys. Rev. D 49, 6410 (1994). 

[35] L. Randall, M. Soljacic and A. H. Guth, Nucl. Phys. B 
472, 377 (1996). 

[36] D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999). 

[37] J. Garcia-Bellido, A. D. Linde and D. Wands, Phys. Rev. 
D 54, 6040 (1996). 

[38] W. R. Gilks, S. Richardson & D. J. Spiegelhalter, Markov 
Chain Monte Carlo in practice, Chapman & Hall: Lon- 
don (1996). 

[39] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 

[40] M. Tegmark et al, Phys. Rev. D 69, 103501 (2004). 

[41] V. Barger, H. S. Lee, and D. Marfatia, Phys. Lett. B 
565, 33 (2003); W. H. Kinney, E. W. Kolb, A. Mel- 
chiorri and A. Riotto, Phys. Rev. D 69, 103516 (2004); 
S. M. Leach and A. R. Liddle, Phys. Rev. D 68, 123508 

(2003) ; S. Tsujikawa and A. R. Liddle, JCAP 0403, 
001 (2004); G. Calcagni and S. Tsujikawa, arXiv:astro- 
ph/0407543. 

[42] Y. S. Piao, S. Tsujikawa and X. Zhang, Class. Quant. 

Grav. 21, 4455 (2004). 
[43] H. Akaike, IEEE Trans. Auto. Control, 19, 716 (1974); 

G. Schwarz, Annals of Statistics, 5, 461 (1978); 

T. T. Takeuchi, Astrophys. Space Sci. 271, 213 (2000); 

A. Nakamichi and M. Morikawa, PhysicaA 341, 215 

(2004) ; A. R. Liddle, arXiv:astro-ph/0401198. 

[44] D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003). 
[45] M. Viel, J. Weller and M. Haehnelt, arXiv:astro- 
ph/0407294. 



